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ABSTRACT 

In this paper we discuss the effect of Coulomb collisions on the temperature profiles of 
the intracluster medium in clusters of galaxies, motivated by recent reports of negative 
temperature gradients in some clusters by Markevitch et al. The timescale for electrons 
and protons to reach temperature equilibrium can exceed a few billion years beyond 
radii of a Mpc, if the intracluster gas is assumed to be at the usual cluster virial 
temperature. If a cluster merger has occurred within that time causing the protons, 
but not the electrons, to be rapidly heated then a small negative temperature gradient 
can result. This gradient is larger in clusters with high temperatures and steep density 
profiles. 

Applying these considerations to the cluster of galaxies A2163, we conclude that, 
more plausibly, the observed gradient is due to a lack of hydrostatic equilibrium fol- 
lowing a merger. 
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1 INTRODUCTION 

Recent results on steep temperature profiles in some clusters 
of galaxies, as measured by ASCA X-ray spectra (Marke- 
vitch 1996 and Markevitch et al. 1996), have raised questions 
on the physical condition of the intracluster medium (ICM). 
At a radius of a few 100 kpc, the ICM has a characteristic 
density n gaB ~ 10~ 3 cm -3 , temperature T gas ~ 10 8 K and a 
heavy element abundance of about 30 per cent of the solar 
value. The density drops at larger radii r approximately as 
r 2 . 

Assuming a polytropic distribution for the gas, the tem- 
perature scales with density as 



(1) 



where the polytropic index 7 ranges between 1 and 5/3, the 
limits corresponding to the gas being isothermal and adi- 
abatic, respectively. Markevitch (1996) found that 7 « 1.9 
and 1.7 for the clusters A2163 and A665, respectively (> 1.7 
and 1.3 and the 90 per cent confidence level). When 7 > 5/3 
the gas is convectively unstable, which tends to produce tur- 
bulent motions in the ICM, resulting in a mixed, homoge- 
neous gas after several sound crossing times. The detection 
of a dramatic drop in the temperature profile may be an in- 
dication that those clusters are observed in a very unusual, 
brief stage in their existence, perhaps having experienced a 
major merger within the previous few billion years. 

An alternative possibility to explain the steep tempera- 
ture decline , suggested by Markevitch (1996 and Markevitch 
et al 1996), is that the electron temperature, T c , which is the 



quantity measured by X-ray observations, is not representa- 
tive of the mean gas temperature, Tgas- This can occur, for 
example if a shock has heated the ICM, with the protons re- 
ceiving most of the energy so that the proton temperature T p 
greatly exceeds the electron temperature T c . Coulomb scat- 
tering then equilibrates the temperatures at a rate (Spitzer 
1962): 



dt 



T n -T„ 



(2) 



In this equation, t eq , the equipartition time via Coulomb 
scattering, is given by 



tc 



n'i/2 



2 m c ripCaT In A ' 



(3) 



where 8 C = kBT e /(m e c 2 ), n p is the proton density, err is the 
Thomson scattering cross section and A is the ratio of largest 
to smallest impact parameters for the collisions (InA ~ 37.8, 
Sarazin 1988; note that a similar term in the corresponding 
dimensionless proton temperature 9 P is negligible for the 
conditions of the ICM). The rate is proportional to the gas 
density and so can be low, and the equilibration time long, 
at the outer parts of a cluster where the density is least. 

Here we examine the equilibration timescale for a typi- 
cal ICM under the conservation of total kinetic energy. We 
find (Section 2) that an implausibly short timescale is re- 
quired for a large temperature gradient to be observed. Both 
the X-ray emission and the equilibration of the ICM depend 
upon Coulomb scattering, so when the rate of emission is 
high the rate of equilibration is also high. In Section 3, we 
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Figure 1. This plot shows the temperature profiles as a function 
of radius r at the time t m of 3 Gyrs. The input parameters are: 
n = 2.5 X 10- 3 cm- 3 ,/3 = 2/3 and r c = 0.3. The thick solid 
line indicate the effect of the e-p Coulomb collisions on the elec- 
tron temperature when the energy density U is conserved and 
assuming that T ga s(?*) is isothermal. The radial decrease of T c 
with respect to the central value is about 0.02, 4 and 14 per cent 
at 1, 2 and 3 Mpc, respectively. The dotted line describes the T p 
profile and the dashed line the projected temperature < T >. 
The latter profile diminishes with respect to the central value by 
1, 9 and 21 per cent, at 1, 2 and 3 Mpc, respectively. 



apply these considerations to the case of A2163 and discuss 
our conclusions in Section 4. 



2 ELECTRON-PROTON COULOMB 
INTERACTIONS 

The accretion of infalling material on the cluster potential 
shock heats the protons to an initial temperature Tp, which 
is representative of their isotropic Maxwellian velocity dis- 
tribution. The protons equilibrate between themselves in a 
time shorter than £ eq by a factor (me/nip) 1 ^ 2 ~ 0.02, so 
that they can be considered to be in equilibrium when the 
heating event occurred a time t m ~ t cq in the past. The elec- 
trons, which are not strongly involved in shock events due to 
their negligible mass, achieve a low temperature T e l <C Tp. 

Since the total kinetic energy density U = §fc(n p T p + 
n e T e ) has to be conserved, we have the following implica- 
tions: (i) the local mean gas temperature, T gas , 



n p T p + n c T c 



T p + 1.21T 
2.21 



(4) 



is constant with time; (ii) the initial Tp is about 2.2 times 
the balanced temperature value, Tg aa ; (iii) T e increases (and 
T p decreases) with time, with the energy exchange between 
protons and electrons driven by the relation 



dT e . , . G?Tp 
_ = _ (np/no) _. 

Rearranging eqn.(2) in the form: 



(5) 



Figure 2. The behavior of the electron temperature as function 
of r/r c , when fi = 2/3 (solid line) and fi = 1 (dashed line). The 
different profile were calculated with different t m (upwards: 1, 3, 
10 Gyrs), assuming constant no and fixing T e (r = 0) at 8 keV. 
The thickest solid line corresponds to fi = 2/3 and t m = 3 Gyrs. 
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Figure 3. The dependence of T c , normalized to the central value 
T c (r = 0), upon different central temperature values (downwards: 
4, 6, 8 - the thickest line - and 10 keV) is shown. All the profiles 
were calculated at t m = 3 Gyrs, assuming fi = 2/3 (solid lines) 
and fi = 1 (dashed lines). 



T p T e 



dT e = 



dt, 



(6) 



and using relation (4), we can solve analytically eqn.(6), 
obtaining 



3 V T Z1 




3/2 



1/2 



cot hiA t m n g 



(7) 
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Figure 4. The dependent variable T c /T gas at first member in 
eqn.(7) is plotted versus the independent variable present in the 
second member. 



Where T e /T gas < 0.5, such as in the outer part of the 
temperature profile given a high central temperature and 
steep density profile, we can further expand the logarithmic 
term in eqn.(6). Then, for T gas =constant, we obtain that 

rp 2/5 

We now parametrise the proton density by a j3— model 
(Cavaliere & Fusco-Femiano 1978), as generally adopted for 
the description of the ICM with central density no and core 
radius r c , 



i p (r) = n \l + {r/r c y 



-3/3/2 



(8) 



and fix the electron density n c = 1.21 x n p . 

Results for reference values of t m = 3 Gyrs, T e (r = 0) 
= 8 keV, n = 2.5 x 10" 3 cm" 3 , r c = 0.3 Mpc and (3 = 2/3, 
are shown in Fig. 1. A significant small temperature gradient 
can remain at this reference time. In Fig. (1), we also show 
the projected temperature, < T >, which is the emissivity- 
weighted T e : 



< T > = 



e(r) T c (r) dr" 



l, 2 



e(r) dr 2 



b 2 ' 



(9) 



where the volume emissivity e is integrated along the line of 
sight at projected radius b. It is more appropriate to com- 
pare with observed values and is shown to decrease more 
significantly than T e (r). 

The dependence of the temperature gradient T e (r) on 
tmt P an d the central value of T e , T c (r = 0), is shown in 
Fig. 2 and 3. Note that the temperature gradient steepen 
when there is (i) a higher energy per particle (i.e. when a 
larger central T e is considered); (ii) a greater value for /3, so 
that the gas density becomes steeper; or, of course, (iii) the 
timescale is shorter. 

In Fig. 4, we show how T e /T sas varies with t m ng as r ga 3/ ' 2 . 
A 10 per cent difference in T^, with respect to T gas , requires 
that n gas < 3.09 x lO -4 ^ 1 ^ 2 cm , where t m is in Gyrs 
and T gas in units of 10 keV. We note that the outer densities 
detectable by the ROSAT PSPC in deep cluster images are 
of order of about 10 _4 cm -3 . 



Figure 5. The temperatures profiles discussed in Sect. 2.1 are 
plotted here, assuming (i) a ratio r c /r200 of 0.15 (cf. table 1 in 
Navarro et al. 1996) for the dotted profile from the N-body simu- 
lations, and (ii) a j3- model with (3 = 2/3 for the adiabatic profile 
(dashed line). The shadowy region indicates where the gradients 
cause the ICM to be convectively unstable. 



2.1 The presence of a gradient in the gas 
temperature 



Above, we have considered the case of T gas (r) = constant. 

However, recent N-body simulations of clusters of galax- 
ies (Navarro et al. 1996, Evrard et al. 1996) have shown that 
T gas is almost constant in the range 0.1-0.4 rjrioo (where 
r2oo indicates the radius at which the cluster mean density 
is 200 times the critical density value) and slightly decreases 
as r- -1//2 towards the outer virialized part. Such a gradient is 
not unstable to convective mixing, from the adiabatic con- 
dition [cf. eqn.(l) and Fig. 5]: 



dT e: 



d 



( 2 /A 

\n s L I 



dr ' dr ' ^ 

We can use eqn.(7) (obtained under the condition of 
T gas constant with time) to investigate the more extreme 
electron temperature gradient when the gas is adiabatic. We 
obtain a behavior of T c almost coincident with the T gas pro- 
file shown in Fig. 5, for every t m considered. 



3 THE CASE OF A2163 

A long Coulomb equilibration timescale between electrons 
and protons has been suggested by Markevitch et al. as one 
explanation for the steep negative temperature profile ob- 
served in ASCA observations of A2163. We now try to test 
this possibility by comparing our calculated profile with the 
deprojected temperature in Abell 2163 (Markevitch et al. 
1996; Elbaz et al. 1995; Fig. 6). In particular, we find that it 
is impossible to reconcile the theoretical temperature gradi- 
ents with the observed one in the outer parts of the cluster, 
for any reasonable i m (even for 1 Gyr the disagreement is 
about of 5.4 a), if the final state for the ICM is isothermal. 
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Figure 6. The temperature profile of A2163 from the analysis by 
Markevitch et al. (1996). The model profiles are calculated using 
as input the values in Elbaz et al. (1995): no = 6.65 X 10 — 3 cm -3 , 
j3 = 0.62. T B (r = 0) is considered both equal to 12.2 (solid lines; 
Markevitch et al.) and 14.6 (dotted lines; Elbaz et al.) keV. The 
results from an assumed isothermal profile clearly fail at the outer 
radius. The two adiabatic profiles (dashed lines) are calculated at 
r/r c equal to and 4.5 via eqn.(l) and with 7 = 5/3. From the 
discussion in Sect. 2.1, they do not depend significantly on t m . 



Only if we make the ad hoc assumption that the ICM is adia- 
batic beyond about one Mpc can we obtain some agreement 
(cf. dashed lines in Fig. 6). 

We now consider whether the gas can be considered to 
be in hydrostatic equilibrium since the last plausible merger 
(a few billion years). The time required for a sound wave to 
cross a cluster is obtained by the relation w| = dP gas /dp gas , 
where P gas and p gas are the ICM pressure and density, re- 
spectively. Solving this relation with respect to time, we get 
ts oc r Tgas' 5 , that does not depend significantly on any tem- 
perature (density) gradient. Using the values in Fig. 6 for 
A2163, we calculate that t s is 0.6±0.6, 2.1±0.8 and 7.7±2.9 
Gyrs at 0.45, 1.43 and 2.99 Mpc, respectively (the errors are 
la deviation obtained by propagation on r and T gas uncer- 
tainties). A critical condition for the hydrostatic assumption 
appears above 1.4 Mpc (~ 5r c ), where the steep tempera- 
ture profile requires that a merger occurred less than about 
3 Gyrs ago, if the steepness is due to the Coulomb pro- 
cess. Then, the cluster appears unrelaxed in its outer parts, 
with the subsequent presence of gas bulk motions, whose 
pressure is not accounted for by the measured ICM temper- 
ature in the last spatial bin of Fig. 6 (see also the discussion 
in Markevitch et al. 1996). We conclude that in regions of 
(r, t m ) space where the electron temperature is significantly 
below the proton temperature then it is likely that hydro- 
static equilibrium has broken down. 

Another timescale relevant to the thermalization of the 
ICM is the ionization equilibrium time of the abundant ele- 
ments, ti n- The emission lines of the highly ionized Fe atoms 
are the strongest features detected in numerous X-ray spec- 
tra of clusters of galaxies. Using the fitted parameters in 
Shull & Van Steenberg (1982) for the collisional ionization 
rate coefficients Cf c of the hydrogenic FeXXVI ions, we get 



iion = (Cf c n c ) _1 ~ 0.1 — 10 Gyrs in the range < r/r c < 10 
for the density profiles considered above (i.e. /3 = 2/3 and 
1, T e = 8 keV). In the case of A2163, ii on is of the order of 
a few tenth of Gyr, and less than 0.4 Gyr, which is shorter 
than any other timescale discussed here. 

Finally, we note that the thermal conduction timescale 
on which temperature gradients are erased are much longer 
than the other timescales considered here by more than one 
order of magnitude. 



4 CONCLUSIONS 

In this paper, we have shown that electron-proton Coulomb 
interactions in ICM become inefficient in reaching the 
equipartition, steepen the electron temperature gradient sig- 
nificantly, only if (i) the energy per particle is high, (ii) the 
gas density profile is steep, and (iii) the time elapsed since 
the last merger, or proton heating event, is very short, i.e. 
n gas < 3.09 x 10" 4 (f m /lGyr)" 1 (r gas /10keV) 3/2 cm- 3 . Local 
conservation of energy means that the proton temperature 
drops as the electron temperature rises. Regions of clusters 
where a large disequilibrium occurs are likely to both be out 
of hydrostatic equilibrium and to have low X-ray emission. 

Similar conclusions are also presented by Fox & Loeb 
(1997), who consider Coulomb interactions in a plasma ac- 
creting on a cluster in the framework of the spherical self- 
similar model. 

When applied to A2163, for which Markevitch et al. 
(1996) report a steep drop in electron temperature in the 
outer parts of the cluster, we find that we cannot reproduce 
the profile from gas in hydrostatic equilibrium, without re- 
quiring the mean gas temperature to drop sharply beyond 
about one Mpc. 
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